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ABSTRACT 

Focusing on the well motivated aperture mass statistics Map, we study the possi- 
bility of constraining cosmological parameters using future space based SNAP class 
weak lensing missions. Using completely analytical results we construct the covariance 
matrix for estimators based on two-point and three-point statistics. Our approach in- 
corporates an accurate modelling of higher-order statistics to describe cosmic variance 
as well as various sources of discrete noise at small angular scales. These results are 
then fed into a Fisher matrix based analysis to study cosmological parameter degen- 
eracies. Joint and independent analysis, with or without redshift binning, for various 
parameter combinations are presented. An analytical modelling of the covariance ma- 
trix opens up the possibility of testing various approximations which are often used in 
derivations of semi-analytical results. These include how inclusion of full non-Gaussian 
terms in covariance matrix affects parameter estimation. Inclusion of three-point in- 
formation and how such information can enhance the accuracy with which certain 
parameters can be estimated is also studied in detail. It is shown that broad corre- 
lation structure among various angular scales in such circumstances mean reduction 
in number of available angular scales which carry completely independent informa- 
tion. On the other hand, the effect of theoretical inaccuracies, in modelling either the 
power-spectrum or bi-spectrum evolution, onto the measure of cosmological param- 
eters from weak lensing surveys is also considered. Several cosmological parameters, 
fim, CTg, spectral index n^, running of spectral index as and equation of state of the 
dark energy Wdc are included in the analysis. 

Key Vifords: Cosmology: theory - gravitational lensing - large-scale structure of 
Universe - Methods: analytical, statistical, numerical 



1 INTRODUCTION 

Weak lensing surveys are expected to make significant con- 
tribution in constraining cosmology by studying the statis- 
tics of dark matter distribution in tlie nearby universe in 
an unbiased way. For instance, Contaldi et al.(2003) used 
the Red Cluster Sequence (RCS) to show that the and 
as degeneracy directions are nearly orthogonal making them 
particularly suitable for combined analysis (van Waerbeke et 
al. 2002). Ishak et al. (2003) argued that joint CMB-cosmic 
shear surveys provide an optimal data set for constraining 
the amplitude and running of spectral index which helps to 
probe various inflationary models. Tereno et. al.(2004) have 
studied cosmological forecasts for joint CMB and weak lens- 
ing survey data. Clearly, the potential of weak lensing sur- 
veys (Mellier 1999; Bartelmann & Schneider 2001; Refregier 
2003; van Waerbeke & Mellier 2003) as a cosmological probe 
is now well established (Contaldi et al. 2003; Hu & Tegmark 



1999). In last few years there have been many studies which 
have directly detected cosmological shear in random patches 
of the sky (Brown et al. 2003; Bacon et al. 2002; Bacon, 
Refregier & Ellis 2000; Hamana et al. 2003; Hammerle et 
al. 2002; Hoekstra et al. 2002a; Hoekstra, Yee & Gladders 
2002a; Jarvis et al. 2002; Kaiser, Wilson & Luppino 2000; 
Maoh et al. 2001; Refregier, Rhodes, & Groth 2002; Rhodes, 
Refregier & Groth 2001; van Waerbeke et al. 2000; van 
Waerbeke et al. 2001a; van Waerbeke et al. 2002; Wittman 
et al. 2000). While early studies were primarily concerned 
with detection of weak lensing signal, present generation of 
weak lensing studies are putting constraints on cosmologi- 
cal parameters such as matter density parameter ilm and 
power spectrum normalisation as- Not only these studies 
can put remarkable constraints on cosmological parameters, 
they can help to break up parameter degeneracies when used 
along with other cosmological surveys such as CMB obser- 
vations as we discussed above. 
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Inspired by the success of these surveys, there are many 
other ongoing, planned and proposed weak lensing surveys 
which are currently in progress such as the Deep Lens Sur- 
vey (Wittman et al. 2002), the NOAO Deep survey (Melher 
et al. 2001) ^, the Canada- France-Hawaii Telescope Legacy 
Survey (Mellier et al. 2001) ^, the Panoramaic Survey Tele- 
scope and Rapid Response System ^ , the Supernova Accel- 
eration Probe " (Rhode et. al. 2003; Massey et al. 2003) and 
the large Synoptic Survey Telescope ^ (Tyson et. al. 2002). 
Independently or jointly these observations can contribute 
to increase the accuracy with which various cosmological pa- 
rameters can be determined from other cosmological obser- 
vation such as CMB. While first generations of cosmic shear 
surveys have demonstrated (van Waerbeke & Mellier 2003; 
Refregier 2003) the feasibility of weak lensing studies in con- 
straining the dark matter power spectrum parametrised by 
(Ts, f^m and shape parameter V (Munshi & Valageas 2005 
contains a recent compilations of present ground based sur- 
vey results), the future surveys will be able to probe much 
larger scales and therefore it will be possible to study the 
linear regime directly to put more stringent bounds on cos- 
mological parameters such as the equation of state of dark 
energy and its time variations. 

Various groups have carried out detailed analysis of 
cosmological constraints that are achievable from ground 
and space based weak lensing surveys (Beruardeau et al. 
1997, Jain & Seljak 1997, Kaiser 1998). Early studies fo- 
cused mainly on a limited number of cosmological parame- 
ters whereas later on it was realised that weak lensing can 
play am important role in constraining the dark energy equa- 
tion of state (Hu 1999, 2002ab; Hueter 2002; Abazajjian & 
Dodelson 2003; Heavens 2003; Refrefier et al. 2003; Benabed 
& van Waerbeke 2003; Jain & Taylor 2003 and references 
therein). Hu & Tegmark (2004) studied a general multi- 
parameter Fisher matrix analysis based on power spectrum 
to check prospects of weak lensing for parameter estimation 
independently or jointly with external data set such as CMB. 
Takada & Jain (2005) have presented a study based on joint 
analysis of power spectrum and bi-spectrum and emphasised 
the role of tomography in in weak lensing parameter estima- 
tion. More recently Kilbinger & Schneider (2005) have stud- 
ied joint constraints from second- and third- order moments 
based on aperture mass Map statistics to study cosmological 
parameter estimation efHciency. Although similar in moti- 
vation we use a different analytical approach to check the 
efficiency of future space based weak lensing surveys to con- 
strain cosmology. We combine measurements from second 
(Mfp) and third order (Mfp) statistics to constrain cosmo- 
logical parameters such as Jim, fs, as, ris and Wde- The effect 
of estimating the median redshif zq of the source redshift dis- 
tribution directly from the data is also included. Impact of 
inaccurate modelling of power-spectrum and bi-spectrum of 
the underlying mass distribution is also analysed. 

The reasons for choosing Map over other two-point 
statistics, as pointed out by Kilbinger & Schneider (2005), 
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are several. Unlike other two-point statistics Map by con- 
struction is a scalar object and unlike many other three- 
point objects, constructed from shear three-point correla- 
tion functions, the third order moment of Map has a non- 
vanishing ensemble average and hence can directly probe 
the underlying bi-spectrum of the matter distribution. Be- 
sides, Map and its extensions has the inherent characteristic 
of separating gravity induced E mode from B mode (see e.g. 
Crittenden et al. 2001, for more detailed discussion on this 
issue) generated mainly by residual systematics which can 
affect both second- and third-order estimators. Moreover, it 
is to be noted that the integral measure (Mfp) is also a lo- 
cal measure of the bi-spectrum and does not suffer from the 
oscillatory nature of other filters. 

This paper is organised as follows. In section 2 we in- 
troduce the second- and third- order estimators based on 
Map statistics. Expressions are presented for their joint co- 
variances. Section 3 describes the cosmological parameters 
and the survey design for proposed SNAP class experiments. 
Next, in sections 4 and 5 we study the influence of non- 
Gaussianities onto the covariance matrices and the signal to 
noise ratios. In section 6 we introduce the Fisher matrix for- 
malism in our context and use it to study estimation error 
for weak lensing observable. Finally the section 5 is left for 
discussion of our results and future prospects. Appendix A 
is devoted to comparing the performance of ground based 
surveys and SNAP class space based surveys. It also takes a 
detailed look at various issues of weak lensing survey design. 

2 APERTURE-MASS STATISTICS 
2.1 Weak-lensing effects 

Weak gravitational lensing effects probe the matter distribu- 
tion integrated along the line of sight to distant sources like 
galaxies (e.g., Beruardeau et al. 1997; Kaiser 1998). There- 
fore, they can be used to obtain information on the 2-D 
projected density field k (also called the convergence field) 
on the sky. In practice, one needs to talce into account the 
redshift distribution n{zs) of the sources. Indeed, in order 
to have good statistics one needs to average weak lensing 
effects over many sources which entails rather broad source 
distributions. On the other hand, it can be convenient to use 
a filtered version of the convergence n, such as the aperture- 
mass Map. More specifically, the latter is obtained from n 
by using a compensated filter. Then, one can also express 
Map as a function of the tangential shear 7t (Kaiser 1994; 
Schneider 1996) which can be directly estimated from the 
observed ellipticities of distant galaxies. This is the property 
which makes the aperture-mass a useful quantity. In addi- 
tion, since it involves a compensated filter the contribution 
from long wavelengths to the signal is damped as compared 
with a top-hat filter (which would yield the mean conver- 
gence) so that A/ap allows one to probe the matter density 
field over a narrow range of wavelengths which can be varied 
through the angular scale 6s of the window Um^p{'&)- Thus, 
the aperture-mass can be written in terms of the fluctuations 
of the density field as: 

Map = j M Um,,0) j^'^^^ Ax w{x) S{x,^&), (1) 
with: 
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dzs n{zs) — 2" 



{l+z).{2) 



Here the redshift z corresponds to the radial distance x a-nd 
2? is the angular distance, i? is the angular direction on the 
sky, 6{x,'D'&) is the matter density contrast and hereafter 
we normalise the mean redshift distribution of the sources 
(e.g. galaxies) to unity: J dzs n{zs) — 1. We defined Zmax as 
the depth of the survey (i.e. n{zs) = for Zs > ^max). Here 
and in the following we use the Born approximation which 
is well-suited to weak-lensing studies: the fluctuations of the 
gravitational potential are computed along the unperturbed 
trajectory of the photon (Kaiser 1992). We also neglect the 
discrete efi'ects due to the finite number of galaxies. They 
can be obtained by taking into account the discrete nature 
of the distribution n(za). This gives corrections of order 1/N 
to higher-order moments of weak-lensing observable, where 
N is the number of galaxies within the circular field of inter- 
est. In practice is much larger than unity (for a circular 
window of radius 1 arcmin we expect TV > 100 for the SNAP 
mission) therefore in this paper we shall work with eq.Q. 
We choose for the filter Um^^ associated with the aperture- 
mass Map the window of finite support used in Schneider 
(1996): 
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where O is a top-hat with obvious notations. The angular 
radius 63 gives the angular scale probed by this smoothed 
observable. As described in Munshi et al. (2004), the cumu- 
lants of Afap can be written in Fourier space as: 



Y[dX^ W(X^) / Y[dk, WM^,{-k±,Vje3 



J = l 



[le"'"'^^' (5(ki)...5(k,)), 



(4) 



We note (..) the average over different realizations of the 
density field, fcy is the component of k parallel to the 
line-of-sight, kj_ is the two-dimensional vector formed by 
the components of k perpendicular to the line-of-sight and 
Wm^p {k±V9s) is the Fourier transform of the window Um^,, '■ 



W^M„p(kxO0s 



(5) 



The Fourier-space expression Q is well suited to mod- 
els which give a simple expression for the correlations 
((5(ki)..5(kp))c, such as the stellar model (Valageas et al. 
2004; Barber et al. 2004) defined by: 



(<5(ki)..5(kp)). 



^£fo(ki + ...+kp) ^J]p(fc,), (6) 



where S2 ~ 1, Sd is the Dirac distribution and P{k) is the 
3-D power-spectrum of the density fluctuations. The coeffi- 
cients S3, S4, .. are closely related (and approximately equal) 
to the skewness, kurtosis, .., of the density field. Eq.Q gen- 
eralises in a straightforward fashion for many-point cumu- 
lants. This allows one to consider for instance the cross- 
correlations between the two statistics Mapi and Map2 asso- 
ciated with two different angular scales dsi,ds2 and/or two 
different redshift distributions ni{zs),n2{zs). 



2.2 Low-order estimators 

The expressions recalled in the previous section describe 
weak lensing effects due to the fluctuations of the mat- 
ter density field. In practice, one actually measures the 
aperture- mass Afap from the distortions of the images of dis- 
tant sources. Thus, in the case of weak lensing {\k\ <^ 1) the 
observed complex ellipticity e — ei + ie2 of a distant galaxy 
is related to the shear 7 = 71-1-172 by: e — 7-fe,, where e» is 
the intrinsic ellipticity of the galaxy. On the other hand, the 
aperture-mass defined in eqs.Q-J^ can also be written as 
a function of the tangential shear 7t as (Kaiser et al. 1994, 
Schneider 1996) as: 



Map= / d^QM^.WjtW 



with: 



(7) 



(8) 



This leads us to define the estimators Mp for low-order mo- 
ments (e.g., Munshi & Coles 2003, Valageas et al. 2004b): 



2\P 



Mr, 



with: 



Mil 



N 

E ' 

(ji,---,jp) 



(N)p = N(N - 1)..{N - p + 1) 



(9) 



(10) 



where A*' is the number of galaxies in the patch of size 7r6's, 
p is the order of the moment and Qj = QM^-p{'&j) where 
•Oj and etj are the direction on the sky and the observed 
tangential ellipticity of the galaxy j. Finally, the sum in 
eq. runs over all sets of p different galaxies among the 
galaxies enclosed in the angular radius 6s , which ensures that 
{Mp) = {Mip) if we neglect the correlation of the intrinsic 
ellipticity of a given galaxy with other galaxies or with weak 
lensing observables. 

These estimators Alp correspond to a single circular 
field of angular radius 9s containing A'^ galaxies. In practice, 
the size of the survey is much larger than 6s and we can 
average over Nc cells on the sky. This yields the estimators 
Aip defined by: 



whence (Mp) = (Mfp), 



(11) 



1 

n = l 
(n) 

where Mp is the estimator Mp for the cell n and we as- 
sumed that these cells are sufficiently well separated so as to 
be uncorrelated. The estimators Mp and Mp provide a mea- 
sure of the moments (Mfp), which can be used to constrain 
cosmological parameters. However, as shown in Valageas et 
al. (2004b), it is better to first consider cumulant-inspired 
estimators Hp. Thus, for the second and third-order cumu- 
lants we define: 

H2 = M2, Ha = M3~3M2Mi, 
and: 



(12) 



1\ c 
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The interest of 7^3 is that its scatter is smaller than for M3, 
see Valageas et al. (2004b). Besides it directly yields the one- 
point cumulants (here we neglected higher-order terms over 

2.3 Covariance matrix 

The estimators Mp and Tip defined in the previous section 
allow one to estimate the cumulants of the aperture-mass 
{Mlp)c for a given angular scale 9s and source redshift dis- 
tribution n{zs). In practice we can vary both the angular 
scale 63 of the filters Um^^ and Qm^^ and the redshift dis- 
tribution n{zs) (for instance by applying a simple binning 
of the galaxies over redshift and selecting different redshift 
bins, which is often referred to as tomography). Thus, if 
we restrict ourselves to second-order moments (Mfp) we ob- 
tain the set of estimators 7^2 (*)i with {Ti2{i)) ~ {M^p{i)) 
as defined in ea. l|18|l . where the index i — l,..,Ni stands 
for both the angular scale 6si and the redshift distribution 
ni{za) (for Ng angular scales and A'^^ redshift bins we have 
Ni = Ne X Nz). Then, the covariance matrix dj associated 
with this data set is: 



a, = {H2ii)H2{j)) - {H2{t))(H2U))- 



(14) 



It measures the cross-correlation between the two estimators 
Ti.2{i) and 7^2 (j). In the following, we shall consider the limit 
where the number of cells Nc on the sky goes to infinity, that 
is the sum in eg. 11311 is performed over Nc cells which cover 
the whole survey area with a uniform coverage and which 
are separated by an angular shift AS which goes to zero. 
Therefore, the discrete sum 11311 tends to the integral: 



■H2ii)= I ^H2{i;d 



(15) 



where A is the survey area and H2{i;a) is the estimator 
H2{i), associated with the angular scale 9si and redshift dis- 
tribution ni{zs), for the cell centred on the direction a on 
the sky. Hereafter we neglect side effects and we consider 
that all estimators cover the same area A, independently of 
the radius Osi. Then, after one angular integration we obtain 
for the covariance Ci, : 



Ci 



A 



a\H2{i),H2(j);a), 



(16) 



where a {H2{i), H2{j); a) is the cross-correlation between 
the two estimators H2(i) and H2{j) separated by the angu- 
lar shift a. From eq.Q the cross-correlation reads (see 
also the expressions given in App.Al of Munshi & Valageas 
2005): 

o\H2{{),H2{j);S) = (Mfp(j)Mfp(j)>,(a) 



-f2[{Map(i)Map(j))c(3) +Q,,(q)]^ 



(17) 



where we recalled explicitly the dependence on the angular 
shift a, and we introduced Qij defined by: 



Qij{a) 



(18) 



where crj — is the scatter of the galaxy intrinsic ellip- 

ticities, rigi is the number surface density of galaxies associ- 
ated with the redshift distribution ni{za), n^ij is the surface 
density of common galaxies to both redshift distributions 



ni{zs) and nj{zs), and Qi is the filter Qm^^ defined in eq.l|HJ 
for the radius 6 si. We can obtain in a similar fashion the 
covariance matrix associated with the data set {7^3(1)}, as- 
sociated with third-order cumulants, as well as the full data 
set {7^2(1), 7^3(1)} where we consider both second-order and 
third-order cumulants. Thus, we have (see also Munshi & 
Valageas 2005): 



a'{H2{t),H3{j)) = (A/4(i)A/4(j)), 



+ 6(Map(i)Mfp(j)>c[(Map(i)Afap(j))c + ( 



(19) 



and: 



a\H3{i),H3{j)) = (Mf 



+9{Mfp(i)Mf 



ap \ 
1-2 



ap(*)M4(j)), 
ap(i))c[(Map(i)Map(j)). 



+ I 



+ 9(M4(i)Map(i)>c(Map(i)Mfp(j)>c 
+ 6[(Map(i)Afap(j))c + Q.,f , 



(20) 



where we did not recall explicitly the dependence on the 
angular shift a. 

It is interesting to estimate the amplitude of various 
contributions to the covariance matrix dj, due to noise, 
Gaussian and non-Gaussian terms. Thus, in addition to the 
full matrix dj described above we also introduce the ma- 
trices Cf/'-''-, C,^'" "- and Cf/-. The matrix C^'" "- only 
includes Gaussian terms which contribute to the cosmic vari- 
ance (i.e. non-Gaussianities and the noise are set to zero). 
This yields for instance for a^{H2{i), H2{j)f'''-''- : 



(21) 



a^{H2{i),H2{j)f-'-'- = 2{Map(i)Map(j)>c. 

We obtain from eg. 1201 a similar expression for 
a^{H3{i),H3{j))'^'''''-, while from eq.^ we see at 
once that a^{H2{i), Hsij))'^''' ''- = 0. Next, the matrix 
^G.s.n. jjjgjy^gg g^jj Gaussian terms which involve the shot 



noise. This yields for a^{H2{i), H2{j)) 



G .s.n. . 



a\H2{i),H2{j)) 



G,s.n. 



4(Map(l)Map(j))cQ»j + 2(3?,. (22) 



We again obtain from eg. 12011 a similar relation for 
(j'^{H-i{i), H-i{j))'^'"-"- , while from eg.llT^ we see at once 
that a2(i?2(j),ff3(j))G.-"- = 0. Note that c'r^''^-^- + C"^'"-"- 
is equal to the matrix dj without non-Gaussian terms. Fi- 
nally, the matrix Cfj"' is equal to the matrix dj where we 
set the noise equal to zero. It includes both Gaussian and 
non- Gaussian terms, and we obtain for {H2{i), Hiij))"-"- : 

a\H2{i),H2{j)r^- = {Mfp(i)Afa'p(j))c + 2{Map(i)Afap(j))?.(23) 

From egs. 11911 - I2OII we obtain similar expressions for 
o-2(i/2(i),i?3(j))"" and G'^{H3{i),H3{j)Y--"-. Note that 
'^^(H2(i),H-j,{j)Y'^' is different from zero. Thus, the com- 
parison between the four matrices dj, Cf^'"^'^' , Cf-'"'^' and 
Cfj"' allows us to evaluate the relative importance of the 
shot noise (associated with the dispersion cr? of the galaxy 
intrinsic ellipticity) and of the cosmic variance (merely due 
to the finite size of the survey), as well as the relative im- 
portance of Gaussian and non-Gaussian terms. 



3 BACKGROUND COSMOLOGY AND 
SURVEY PARAMETERS 

We describe in this section the background cosmology and 
the specific correlation hierarchy for the matter distribu- 
tion which we use in numerical computations. We also focus 
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on the SNAP observational strategy to estimate the accu- 
racy which can be obtained from such weak lensing surveys. 
However, as mentioned before the basic formalism developed 
here remains completely general and specific details studied 
here only serve illustrational purposes. 



3.1 Cosmological parameters 

For the background cosmology we consider a fiducial ACDM 
model with Qm ~ 0.3, Jldc ~ 0.7, wdc = —1, Ho — 70 
km/s/Mpc and as = 0.88. We note f2m the matter density 
(cold dark matter plus baryons) and f^do the dark-energy 
density today at jz = 0. We parameterize its equation of 
state by Wde ~ Pdc/Pdc- For Wdc = — 1 this corresponds to 
a simple cosmological constant. In the following we shall 
investigate the dependence of weak-lensing observable on 
the cosmological parameters flm (keeping ilm + ^dc = 1 for 
a flat universe) and Wdc. We always consider Wdc to be a 
constant independent of time. Then, the Hubble expansion 
rate H{z) reads (e.g., Linder & Jenkins 2003): 



H{z) = ffo Van(l + + f7dc(l + ^)3(i + »dc), (24) 



where a{t) = {l + z)~^ is the cosmological scale factor while 
the dot denotes the derivative with respect to physical time 
(we used f2m-|-f2de = 1). The linear growth factor D{t) obeys 
the differential equation: 



D + 2HD - ^S^mifo (1 + zfD = 0, 



(25) 



which yields for the linear growth rate g{a) relative to a 
critical-density universe g — D/a: 



d^g , 5a.. + (5 - 3mde)ndca"^"'<^'= d<; 



+ 



d In a2 fim -I- ridea-^^do 

(3 - 3u-dc)ndoa"^"'"" 



d In a 



+ - 



-9 = 0- 



(26) 



n,-n + fldca-^'^d^ 

We obtain g{a) by solving numerically the differential 
eg.ra. 

For the matter transfer function T{k) we use the fitting 
formulae provided by Eisenstein & Hu (1999), with f2b — 
0.047, Us = 1 and as = 0. Here is the power-law index of 
the primordial power-spectrum and Us is the running index. 
More precisely, we write the linear matter power-spectrum 
as in Spergel et al. (2003): 

PUk) cx T\k), (27) 

where ko — 0.05 Mpc~^. Thus, the local slope n of the pri- 
mordial spectrum is: 

_ dlnPprim _ , f k \ dn 



dink 



Os In 



(fco) 



dlnfc 



(28) 



Next, we use the fit given by Peacock and Dodds (1996) 
to model the non-linear evolution of the matter power- 
spectrum due to gravitational clustering. In order to in- 
vestigate the sensitivity of weak lensing data on this non- 
linear extrapolation we introduce an additional parameter 
/2. Thus, we modify eq.(21) of Peacock & Dodds (1996) as: 



/NL(a:) 



1 -I- Bf3x + [Ax] 



al3 



1 + {[Ax]''g^Q,^)/[f2Vxy^]y 



1//3 



(29) 





f,=0.85, 1, 1,15 



1 10 100 
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Figure 1. Left panel: The power A'^(k)/k at redshift z = 1 
as a function of the comoving wavenumber k. The three upper 
curves correspond to /2 = 0.85 (lower dashed curve), /2 = 1 
(middle solid line, this corresponds to Peacock & Dodds 1996) 
and /2 = 1.15 (upper dotted curve). This corresponds roughly 
to a 10% variation of the matter power-spectrum in the highly 
non-linear regime. The lower dotted curve shows the linear power 
A|^(fc)/fc. Right panel: The variance {M^p) of the aperture-mass 
as a function of smoothing angle 8s for the SNAP survey. We 
again display the cases /2 = 0.85 (lower dashed curve), /2 = 1 
(middle solid line) and /2 = 1.15 (upper dotted curve). 



This is identical to the formula from Peacock & Dodds 
(1996) for /2 = 1 (see their paper for the meaning of various 
parameters). In the linear regime we merely have /NL(a;) — x 
while in the non-linear regime we have fNhix) oc fix?'^'^. 
Therefore, the parameter /2 describes the amplitude of 
the non-linear part of the matter power-spectrum as com- 
pared with Peacock & Dodds (1996) (which corresponds to 
h = 1). 

For illustration purposes, we show in Fig. the non- 
linear power A^(fe)/fc (left panel, as a function of comoving 
wavenumber k) at redshift 2 = 1 and the variance (Mfp) 
(right panel, as a function of smoothing angle ds) for the 
SNAP survey described in section 13.21 below with no red- 
shift binning. We display the curves obtained for the cases 
/2 ~ 0.85, 1, 1.15 from bottom to top. This corresponds to a 
variation of about ±10% for the power-spectrum P{k) in the 
non- linear regime. Thus, the range /2 — 1 ± 0.15 describes 
roughly the current uncertainty on the non-linear power- 
spectrum in the range of interest. Note that /2 is not merely 
a multiplicative factor to the non-linear power-spectrum as 
the recipe from Peacock & Dodds (1996) also involves a 
rescaling of length scales. Here we defined the 3-D power 
per logarithmic interval A^(fc) = ATvk^ P{k, z). Because of 
the integration along the line of sight the power relevant for 
weak lensing observables is actually A^(fc)/fc, and the typi- 
cal wavenumber associated with the angle Os is k A/(T>6s), 
see Munshi et al. (2004). The lower dotted curve in left panel 
of Fig. Qshows the linear power A|, (fc) /fe. We can check that 
at low wavenumbers and at large angular scales which probe 
the linear regime the dependence on /2 disappears. 

For higher-order correlations of the matter density field 
we use the stellar model © described in Valageas et al. 
(2004a). This is actually identical to the minimal tree-model 
up to third-order moments (Munshi et al. 2004). In partic- 
ular, for the skewness S3 of the 3-D matter density field we 
interpolate between the quasi-linear limit S'^^ and the non- 
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Figure 2. The skcwness Sg^"^ of the aperture-mass as as a func- 
tion of tlie smoothing angle 9s for the SNAP survey. We display 
the cases /a = 0.5 (lower dashed curve), /a = 1 (middle solid 
line) and /s = 1.5 (upper dotted curve) which correspond to a 
50% variation of the skewness in the highly non-linear regime. 

linear prediction Sf^^'^ of KEPT (Scoccimarro et al. 1998). 
In order to investigate the sensitivity to the non-linear fit 
used for S3 we introduce a parameter /s and we define the 
skewness Sf^ reached in the highly non-linear regime as: 



Then, on linear scale we use S3 — S'^^ while on highly non- 
linear scales we use S3 = S^^ , see Munshi et al. (2004) for 
details. 

We show in Fig. H the skewness sf^^ = {M%)/{Mlpf 
as a function of the smoothing angle 9s for the cases /a = 
0.5, 1, 1.5, for the SNAP survey described in section be- 
low with no redshift binning. Thus we now consider a varia- 
tion of ±50% on /s and S3 since the skewness is not known 
to better than 50% in the highly non-linear regime (the ac- 
tual current uncertainty may actually be even larger). Of 
course, since the aperture-mass is linear over the density 
fluctuations within our weak-lensing approximation Q this 
yields a 50% variation for st,''^^ at small angular scales. We 
can check again that at large angular scales which probe the 
quasi-linear regime the dependence on /a disappears. 

Thus, the parameters /2 and /a allow us to evaluate 
the sensitivity of weak lensing results onto the amplitude 
of the two-point and three-point matter correlations in the 
non-linear regime, where they are not known to up to high 
accuracy yet. Unless otherwise stated, we use /2 = 1 and 
/a = 1 which correspond to the fits obtained from Peacock 
& Dodds (1996) and from Scoccimarro et al. (1998). 

3.2 Survey parameters 

Hereafter, we adopt the characteristics of the SNAP mission 
as given in Refregier et al.(2004). More precisely, we consider 
the "Wide" survey where the redshift distribution of galaxies 
is given by: 

n{zs) cc Zs e"'^'^^"-* and zq = 1.13, Zmax = 3. (31) 

The variance in shear due to intrinsic ellipticities and mea- 
surement errors is a, = {|e, p)^/^ = 0.31. The survey cov- 
ers an area A — 300 deg^ and the surface density of us- 
able galaxies is Ug — 100 arcmin"^. In order to extract 
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Figure 3. The logarithmic derivative d ln{Mlp) c / d In X of low 
order cumulants with respect to cosmological parameters X. Up- 
per left panel: The dependence of the variance (M^p) on flm (solid 
line), erg (dashed line), tOde (dotted line), Ug (dot dashed line) and 
as (dotted line), as a function of angular scale 6s with no redshift 
binning. For as we plot d \n(M^p) c / das since Os = in our fidu- 
cial model. Upper right panel: Same as upper left panel for the 
third-order cumulant (M^p)c- Lower left panel: Same as upper 
left panel but for two redshift bins. The low redshift bin is shown 
by the same line styles as in upper left panel whereas the high 
redshift bin is shown by the long dashed line. Lower right panel: 
Same as lower left panel for the third-order cumulant {M^p)c. 

some information from the redshift dependence of weak lens- 
ing effects we also divide the "Wide" SNAP survey into 
two redshift bins: < Zs < z, and z, < Zs < Zmax- We 
choose z, — 1.23, which corresponds roughly to the sepa- 
ration provided by the SNAP filters and which splits the 
"Wide" SNAP survey into two samples with the same num- 
ber of galaxies (hence rig = 50 arcmin"^). Note that one 
cannot use too many redshift bins at it decreases the num- 
ber of source galaxies associated with each subsample (for 
the aperture mass we could still obtain good results with 
three bins but we shall restrict ourselves to two redshift bins 
in this paper). The redshift bins that we use are similar to 
those which were used by Refregier et al.(2003) using pho- 
tometric redshifts, except that they have a sharp cutoff and 
non-overlapping source distributions. Note that using over- 
lapping source distributions (over redshift) would increase 
the cross-correlations. 



4 DEPENDENCE OF LOW-ORDER 
ESTIMATORS ON COSMOLOGY 

We show in Fig. 03 the logarithmic derivative 
91n(M|'p)c/91nX of cumulants of order p = 2 (left 
panels) and p = 3 (right panels) with respect to cosmo- 
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logical parameters X, with no redshift binning (upper 
panels) and with two redshift bins (lower panels). For the 
cosmological parameter as we merely plot d\n{Mip)c/da3 
since Os = in our fiducial model. The rather large values 
obtained for these derivatives shows that weak lensing 
effects could potentially constrain cosmological parameters 
up to a good accuracy. By contrast, the small derivative 
with respect to Wdc shows that the equation of state of 
the dark energy component cannot be measured up to a 
similar accuracy. As expected, since weak lensing effects 
are proportional to matter density fluctuations, eq.Q, the 
derivatives with respect to Qui and as are positive over 
most of the angular range 0.1' ~ 1000'. In particular, we can 
check that (Mfp)c ~ (t| as expected from eq.Q. On the 
other hand, we can check that d(Mlp)c/dns > at small 
scales which probe high wavenumbers k and it crosses zero 
at about 20' which corresponds to the normalisation of the 
power-spectrum. Since these derivatives show different vari- 
ations with the angular scale 8s for different cosmological 
parameters, one should be able to constrain simultaneously 
these parameters by using several angular scales. Note that 
by looking over such a large range of angular scales one 
can discriminate the behaviours of different cosmological 
parameters which helps to remove degeneracies. On the 
other hand, one can use the additional information pro- 
vided by higher-order cumulants, such as the third-order 
cumulant (whence the skewness) shown in the right panels. 
However, we can see that for most parameters {M^p)c 
behaves roughly as (M^p)^, except for Qm where there 
is a residual dependence which can be used to measure 
for instance both Sim and as (see also Bernardeau et al. 
1997, as well as Kilbinger & Schneider 2005). Alternatively, 
one can split the survey into two redshift bins and take 
advantage of the different dependence on cosmology of weak 
lensing effects associated with each bin. This is shown in 
both lower panels. Although the curves are quite similar we 
shall check in sect.|S|that using such a redshift binning does 
indeed improve the constraints on cosmological parameters. 



5 COVARIANCE MATRICES OF LOW-ORDER 
ESTIMATORS 

In order to use the estimators Hp to measure cosmologi- 
cal parameters, through their dependence on cosmology dis- 
played in Fig. 13 we need the covariance matrices dj intro- 
duced in sect. 12.31 The latter are necessary to obtain the rel- 
evant error bars through a Fisher matrix analysis fsect. l^TTll 
or a likelihood function f sect. 1^31 . 

5.1 Impact of noise and non-Gaussianities 

We first show in Fig.|l]the covariance matrix dj of the esti- 
mator 7i2{0s) of the variance (Mfp), see egs. 1131 - 11411 . along 
the diagonal i ~ j with no redshift binning {9si — dsj)- 
The full solid line shows the full covariance C while other 
line styles correspond to the specific contributions (j^.s.n.^ 
C^'" "- and d "-, defined in eqs.^-^. As expected, we 
check that at small angular scales 8s < 1' the covariance C 
is dominated by the shot-noise due to the galaxy intrinsic 
ellipticity dispersion, whereas at large scales da > 2' the co- 
variance is dominated by the cosmic variance (i.e. the error 




0.1 1 10 100 1000 

6^ [ arcmiri ] 

Figure 4. Covariance of the estimator 'H2{6s) of the variance 
(Mfp), eqs.C3-OJ- The solid line shows the full covariance 
Cii = cr'^ {'H2{d si) ,'H2{d si)) as a function of smoothing angular 
scale 9 si, with no redshift binning. The dotted line displays the 
contribution C'^'" "' from the galaxy intrinsic ellipticity disper- 
sion to the Gaussian part of C. The dashed line C'G.c.u. ^j^g 
cosmic variance contribution to the Gaussian part of C while 
(jc.v. jg ^j^g cosmic variance contribution to the full matrix C. 




6^ [ arcmin ] 

Figure 5. Covariance of the estimator Ti.3{6s) of the third-order 
cumulant {M^p)c, eg. lliji . The various line styles show different 
contributions to the full covariance Ca = {'H3{dsi),'Ha{9si)) 
(solid line) as in Fig. El 

bar due to the finite size of the survey). We can also check 
that the shot noise contribution C grows as l/e? at 
smaller scales, which is proportional to the number of dis- 
tinct patches of radius 9^ in the survey. The slow rise of the 
cosmic variance term C"'^' yields a broad plateau for the full 
covariance, from 1' up to 1000' , while non-Gaussianities be- 
come negligible above 30'. We can note that even at smallest 
scales non-Gaussian terms are only twice larger than Gaus- 
sian terms so that the Gaussian part C'^ = (7C;,s.n. _^(-,g,c.«. 
of the covariance matrix is always a good approximation 
along the diagonal {8si = 8s j) for two-point estimators such 
as 7^2- 

Next, we display in Fig. |3 the covariance matrix Ca of 
the estimator Tl:i{6s) of the third-order cumulant (Map)c, 
again along the diagonal i = j with no redshift binning 



8 Munshi & Valageas 





C.,(M=,) \_ 






\c 












■ CO.,. 









10 100 1000 

a rem in 1 



1 

r c 


1 \ 


1 i.|-Mnll 


1 1 



10 100 1000 

arcmin | 



: _ C 


C.,(M=,) - 


1 


















? 0.1 















0.1 1 10 100 

9^, [ arcmin ] 



e.-iooo' 



1 10 100 

Sj, I arcmin | 



Figure 6. Covariance dj of the estimator H2{dB) of tlie variance 
(Mfp), egs. 1131 - 1141 . for different angular scales as a function of 
Ogi at fixed dgj , with no redshift binning. The various line styles 
show different contributions to the full covariance as in Fig.|l] The 
four panels correspond to 9sj = 0.1' (upper left panel), 8gj = 2' 



(upper right panel), B^j 
(lower right panel). 



: 50' (lower left panel) and 9^ 



1000' 



{9 si = dsj), see eg. 12011 . We find again that the shot-noise 
due to the galaxy intrinsic ellipticities dominates below 2'. 
However, it now grows as l/9i rather than 1/6^ at smallest 
scales (this is due to the term Qfj in eq. CT). Above 1' the 
covariance is dominated by the cosmic variance but in con- 
trast with Fig. 01 we now find that non-Gaussian terms are 
larger than Gaussian terms by an order of magnitude around 
9a ~ 3' (for the cosmic variance part). This is not really sur- 
prising since Tls itself is an estimator of non-Gaussianities. 
However, above 60' non-Gaussiannities become negligible 
again as we probe the quasi-linear regime. 

We show in Fig. |S] the covariance dj of the estima- 
tor Ti.2{9s) of the variance (Mfp) for different angular scales 
{9si,9sj} with no redshift binning (i.e. we probe dj out 
of the diagonal whereas Fig. 21 was restricted to the diag- 
onal i = j). The line styles are as in Fig. 0] The wiggles 
in the lower panels for C^''^'"' correspond to a change of 
sign (hence we actually plot ^(j'3,s.n.^-^ because at large and 
different angular scales the first term in eg. (1221 can make 
QG,s.n. jiegative.When both angular scales are small the co- 
variance is dominated by the intrinsic ellipticity noise (both 
upper panels with < 2') but as soon as one scale is larger 
than 10' the covariance is dominated by the cosmic variance 
(both upper panels with 9si > 2' and both lower panels for 
all 9si). Then, except at large scales where 9si ~ 9sj the 
cosmic variance is dominated by the non-Gaussian terms 
which can be several orders of magnitude larger than the 
Gaussian contribution. Note that this is quite different from 
Fig. 21 which showed that non-Gaussian contributions were 
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Figure 9. Inverse of the signal to noise ratios of the esti- 
mators 7i2 and 'Hs of the second-order and third-order cu- 
mulants of the aperture-mass. The solid lines show the ratios 
a{Hp{i), Hp(i))/ Hp{i) as a function of the angular scale B^i for 
p = 2 (curve labeled H2) and p = 3 (curve labelled H3) with no 
redshift binning. The covariance = Cu also corresponds to the 
solid line in Figs. IUSI The dotted lines show the ratios obtained 
when we only include Gaussian terms in the covariance Cu while 
the dashed lines correspond to the uncertainty associated with a 
15% variation of /2 with respect to 7^2 and a 50% variation of /a 
with respect to Ha . The dot-dashed line which follows closely the 
curve obtained for H.^ shows the noise/signal ratio when we use 
the estimator fea. llll ') instead of H3. 



not very important along the diagonal. Moreover, the non- 
Gaussian contribution to dj shows a broad plateau as we go 
farther from the diagonal. This feature, which makes the co- 
variance matrix very broad about the diagonal and reflects 
a strong correlation between various angular scales, leads 
to difficulties for the estimation of cosmological parameters 
since the covariance matrix dj cannot be easily inverted. 

The full covariance matrix along with its various compo- 
nents are shown in Fig.|7|for 0.2. In agreement with Figs.^]- 
|S]the upper right panel which displays C^'"-"- shows that 
the shot noise increases fastly at smaller angular scales. Its 
contribution to the covariance is restricted to the diagonal 
9ai = 9s2 above 10' while it mostly depends on max(0si, Ss2) 
when both scales are below 10'. On the other hand, the 
Gaussian contribution C'^''^'^' to the cosmic variance (lower 
right panel) is always restricted close to the diagonal and in- 
creases at large scales. The upper left panel which displays 
C"^"' clearly shows that non-Gaussianities bring a signifi- 
cant broadening to the covariance matrix which is no longer 
diagonal-dominated. Finally, the lower left panel which dis- 
plays the full covariance matrix C refiects these various fac- 
tors and shows a very broad shape with a steep increase 
at small angular scales due to the shot-noise. Fig. m which 
shows the covariance matrix for Ti.^ exhibits a similar be- 
haviour albeit with a stronger impact of non-Gaussianities. 



5.2 Signal to noise ratios 

Next, we plot in Fig. |n| the inverse N/S of the signal to 
noise ratios of the estimators Ti.2 and Tis of the second- 
order and third-order cumulants of the aperture-mass. We 
show the ratios a{Hp{i) , Hp{i)) / Hp{i) without redshift bin- 
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Figure 7. Covariance matrix for the estimator W2 for different angular scales. We plot contours of constant logC Lower left panel: 
the full covariance C. Upper left panel: the cosmic variance contribution C"^"' (i.e. the shot noise due to the galaxy intrinsic ellipticity 
dispersion is set to zero). Lower right panel: the Gaussian contribution C'C.c.u. ^j^g cosmic variance. Upper right panel: the shot noise 
contribution C'G.s.n. ^^i^ Gaussian part of C. 

H3H3-1 bin 




Figure 8. Same as previous figure but for for the estimator H3. See text for discussion. 
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ning (solid lines). They grow at small scales because of the 
shot noise due to the galaxy intrinsic ellipticity dispersion 
and at large scales because of cosmic variance. We also dis- 
play the noise to signal ratios N/S obtained when we only 
include the Gaussian terms C'^ in the covariance Cu (dotted 
lines). This corresponds to = c"'" "- + C'^'^"-, that is 
the sum of dotted lines and dashed lines shown in Figs. I4I5I 
In agreement with sect. 15.1 1 we find that non-Gaussian terms 
only make a difference in the intermediate range 1' — 100' and 
that this effect is quite modest for the estimator Ti.2 of the 
variance (Mfp). For Tis the noise ratio can be increased by 
a factor 3 around 3'. As discussed in sect . 15. ll non-Gaussian 
terms are mainly important for cross-correlations between 
different angular scales. 

It is interesting to note that the angular scale at which 
the highest signal/noise ratio is achieved shifts to smaller an- 
gular scales when we include the contribution to cosmic vari- 
ance from non-Gaussianities. Thus ignoring non-Gaussian 
terms in cosmic variance not only gives a wrong impression 
of higher signal to noise ratio, it also gives a wrong estimate 
of the angular scale where this is achieved. More detailed 
analysis of how these plots change with variations in survey 
charecteristics are presented in the appendix, see fig. lAll 
and fig.m The very high S/N achieved by surveys such as 
SNAP for Ti.3 clearly indicates the potential of weak lensing 
surveys in studying even higher order non-Gaussianity such 
as 7^4 which estimates the kurtosis of underlying mass distri- 
bution. The dot-dashed line which follows closely the curve 
obtained for Tis shows the noise/signal ratio when we use the 
estimator M3 (ea. illl ) instead of Tis. Thus, in agreement 
with Valageas et al.(2005) we find that for third-order esti- 
mators there is not much gain to be obtained by switching 
from Ms to Hs- However, for higher-order non-Gaussianities 
like Ti.4 the improvement can be quite significant (see Figs. 4 
and 5 of Valageas et al. 2005). 

Finally, the dashed lines in Fig.|U]show the noise associ- 
ated with the current uncertainty on the non-linear regime 
of gravitational clustering, as described by the parameters /2 
and fs. Thus, the curve labeled /2 shows N/S = ATi.2/'H2 
with AH2 = {H2{f2 = 1.15) - H2{f2 = 0.85))/2 associ- 
ated with a 15% variation of /2 as described in Fig. In 
agreement with Fig. Q and Fig. |U] this corresponds roughly 
to a 10% uncertainty on the non-linear power-spectrum. We 
can see that this is actually the dominant source of noise 
over the range 0.2' — 20'. Next, the curve labelled /s shows 
N/S = AHa/H-i with AH3 = (Ha (/a = 1.5) - Hsifs = 
0.5))/2 associated with a 50% variation of fs as described in 
Fig. 121 In agreement with Fig. Inland Fig.|5]this corresponds 
roughly to a 50% uncertainty on the non-linear skewness and 
third-order cumulant. This is the main source of noise over 
the range 0.3' — 20'. Thus, in agreement with some previ- 
ous studies we find that current theoretical uncertainties on 
non-linear gravitational clustering are the limiting factor to 
constrain cosmological parameters from large weak-lensing 
surveys with characteristics similar to the SNAP mission. 



5.3 Cross-correlation of redshift bins 

We show in Fig. ll2l the noise/signal ratio we obtain for each 
redshift bin when we divide the sample into two sub-samples 
separated by the source redshift Za = 1.23. The curves are 
similar to those displayed in Fig. |S| for the full sample. For 
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Figure 12. Inverse of the signal to noise ratios of the estimators 
H2 and Hs as in Fig. I^Jbut for two redshift bins. The solid lines 
correspond to the low redshift sub-sample (zs < 1.23) whereas the 
dashed lines show the results for the high redshift bin (23 > 1.23). 
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Figure 13. Inverse of the signal to noise ratios of the estimators 
'Hii,'H2i and 'H12 associated with the cross-correlation between 
both redshift bins. The solid lines are the full noise/signal ratios 
while the dotted lines are the results which we obtain when we 
only keep Gaussian terms in the covariance matrices. 



clarity we do not show the curves associated with the the- 
oretical uncertainties (/2 and /a) or with the Gaussian ap- 
proximation to the covariance matrices. They also follow the 
behaviour of Fig. El We can see that the signal to noise ra- 
tio is larger for the high redshift bin (note that we plot its 
inverse N/S). This is expected since the longer line of sight 
leads to a larger amplitude of weak lensing effects. How- 
ever at very large angles the signal to noise ratio becomes 
slightly better for the low redshift bin for the estimator Tia 
of non-Gaussianities because of the growth of non-Gaussian 
gravitational clustering. 

Finally, we show in Fig. 1131 the noise/signal ratios we 
obtain for the estimators Tin, 7^21 and 7ii2 which directly 
measure the cross-correlation between both redshift bins. 
Here we defined Hpq in a manner similar to Tip introduced 
in sect. 12.21 such that their mean is: (Tipq) = (A/fp^M^p2)c 
where Mapi (resp. Map2) is the aperture-mass associated 
with the low (resp. high) redshift bin, see Munshi & Valageas 
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Figure 10. Covariance and cross-correlation is displayed for two different redshift bins. Left panel shows the covariance of ^2 for lower 
redshift bin and right panel corresponds to higher redshift bin. Middle panel panel shows the cross-correlation between both redshift 
bins. 
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Figure 11. Same as Fig. ESI but for Hs. 



(2005). Of course we can check that the behaviour of differ- 
ent curves is similar to that obtained in previous figures. In 
particular, the signal to noise ratio is slightly better for 7ii2 
than for 7^21 since the former gives more weight to the high 
redshift bin, in agreement with Fig. 1121 However, we note 
that the inclusion of non-Gaussian terms in the covariance 
matrices now makes a large difference for the third-order 
estimators Ti.2i and 7ii2. Indeed, a Gaussian approximation 
could overestimate the signal to noise ratio by up to a factor 
ten. 

The full covariance matrix is displayed in Fig. 1101 for 
7^2 and Fig. 1111 for Tia. In agreement with the discussion 
above we find that the covariance is larger for the high- 
redshift bin, together with the weak lensing signal itself. 
As seen in Fig n2l the overall effect is still to improve the 
signal to noise ratio for the gigh-2 bin. The shape of the 
covariance matrix within each bin (left and right panels) is 
similar to the behaviour obtained for the full survey with no 
redshift binning (Figs. I7I8II . The cross-correlation between 
both bins (middle panels) shows a similarly broad feature 
due to non-Gaussianities. However, there is no rise at small 
angular scales due to shot noise because there are no com- 
mon galaxies to both bins (hence the shot noise contribution 
vanishes). 



6 ESTIMATION OF COSMOLOGICAL 
PARAMETERS 

6.1 Formalism 

Baye's theorem (see e.g. Jaynes, 2003) provides an interest- 
ing starting point for most parameter estimation studies. 
Assuming a cosmological data vector a the likelihood func- 
tion of the cosmological parameter set can be described 



£(0|a) 



£{&,&)£{&) 
£(a) • 



(32) 



Here £.{0) indicates the prior likelihood function of the pa- 
rameter set Q and £(0|a) denotes the conditional likeli- 
hood function. Normalisation determines £(a) and the prior 
comes from other cosmological observations or as in our case 
it is assumed to be constant as we do not include any other 
observational information. The factor £(a|©) describes the 
distribution function of the observed data vector a for a 
given cosmological parameter 0. Assuming a multi-variate 
Gaussian distribution one can express £(a|©) as 



£(a|0) 



v/(27r)^detC(e) 



exp 



a C a 
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(33) 
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where is the inverse of the covariance matrix C = 

{sJ'a.) — (a"^)(a) of the data vector a being used. Associ- 
ated log-hkehhood statistic is defined as where x'^ = 
a''"C~''"a. The covariance matrix C is a function of under- 
lying cosmological parameters. Its derivatives w.r.t various 
cosmological parameters along with the derivatives of the 
mean fia are computed numerically while constructing the 
Fisher matrix Fai3. 

Assuming a fiducial cosmological model, for a given data 
vector, the estimation error of cosmological parameters asso- 
ciated with an unbiased estimator can be constructed from 
the Fisher matrix formalism (see Tegmark, 1997, Matsubara 
& Szalay 2002, Takada & Jain 2003). The Fisher information 
matrix (Kendall & Stuart 1969) which is related to the the 
inverse covariance matrix for the cosmological parameter, 
assesses how well the data vector can distinguish the fiducial 
model from other models. Thus the error on the parameters 
9a obeys: (A^aAS^) > F~0 . The left hand side is computed 
for a specific values of the parameter vector Oq. The equality 
is obtained for maximum likelihood estimates of the param- 
eters. This inequality - which is also known as Cramer-Rao 
inequality - provides a minimum variance bound for unbi- 
ased estimators and can be used to study estimation error 
and their cross-correlation for various parameter sets for a 
given survey strategy. Analytical expression for the Fisher 
matrix F in terms of the covariance and mean of the data 
can be written as (Tegmark et al. 1997): 



Fo/3 = + ^Tr [(In C)„(ln C)^] 



(34) 



where (In C)a, ~ Ca is the derivative of InC w.r.t the 
parameter Oa, and fj,a denotes the derivative of ^ = (a) 
w.r.t. the parameter a. 

The first term corresponds to the case when only mean 
is being estimated from the data vector whereas the sec- 
ond term is associated error for variance estimations from 
the data vector a. The Fisher matrix is a positive-definite 
matrix. It is dominated by the linear order term related to 
estimation of mean of the data. With other cosmological 
data sets where the mean is fixed, such as mean tempera- 
ture of the CMB, the variance term can be the dominant 
term. Here we include both terms in our analysis. 



F<.,3(eo) 



^91n£(a|e) dlnC{a\&) 



dHnC{a\e) 
dQadQfj 



(35) 



Where (...) represents averaging over all possible data 
realization for fixed cosmological parameter values Q — 
©0- The equations of ellipsoid defined by the equation 
J2a,i3 AOaFaiiAe^ = define regions in a multidimensional 
parameter space which can be interpreted as error bounds 
for an experimental setup. A specific choice for A defines 
a Xa confidence level. Although strictly speaking such an 
interpretation is valid only when the likelihood function is 
Gaussian, for a mildly non-Gaussian case it can still provide 
a valuable idea about errors associated with estimators. 

Therefore, up to a normalisation constant the 
marginalised two-dimensional likelihood function C{da,0i3) 
for two parameters 6^ and 9f3 can be expressed as: 



' exp 



-^(^^-^^"H F-C F- 



a/3 
I3l3 



Here {F~^)ap represents the Q/3-element of the inverse 
of the original higher dimensional Fisher matrix F. The 
marginalised error-ellipses therefore can be directly deduced 
from this expression. The likelihood contours start to deviate 
from their Fisher counterparts when higher-order correction 
terms start to become important at a large distance from 
the minima. 

In general when a partial set of parameters is being 
marginalised over the other parameters the error covariance 
of the remaining parameters is given by a sub-matrix of the 
Full inverse Fisher matrix F~^. Marginalisation in general 
causes a broadening of error ellipses due to reduced level of 
prior information being used. 



F = 



Faa 

F4 B 



Fab 
Fbb 



(37) 



(36) 



It can be shown that the marginalised ua x tia Fisher ma- 
trix F can be expressed as F = Faa — FBsF^^F^g. The 
second term in this expression provides the correction to the 
first term, due to marginalisation. In the absence of cross- 
correlations the correction term vanishes. 

In our study, the data vector a consists of various 
choices of measurements of Ti.2 and Tis at different angular 
scales for different redshift bins a^ — {Ti.2{i\ Za),'Hz{i\ Za))- 
Here we have introduced additional index a for parametris- 
ing the redshift bin Za under consideration. A slightly com- 
pact notation was used in previous sections where i was 
allowed to run over different angular scales and different 
redshift bins. Accordingly the Covariance matrix has a 
block structure with various blocks denoting covariance of 
7^2(1; Za) and 7i2(i; Za) along with their cross-covariances be- 
tween same and different redshift bins. In our analysis we 
have considered various combinations of the estimators 7^2 
and 7^3 independently and jointly and for one or two red- 
shift bins to study errors in parameter estimation. Different 
cases that we have considered can be summarised as below. 

• Various angular scales 6s with no redshift information, 
and only 2-point information from Ti.2(i)- 

• Various angular scales 6s with no redshift information, 
but 2-point Ti.2{i) and 3-point Tizii) information. 

• Various angular scales 6s with redshift information, but 
only 2-point 0.2(1, Za) information. 

• Various angular scales 6^ with redshift information, with 
both 2-point and 3-point Ti.2{i, Za),Ti.3{i, Za) iifformation. 

For each of these choices we analyse various combina- 
tions of cosmological parameters independently or jointly 
with or without a prior information regarding median source 
redshift zq. Without prior redshift information the joint co- 
variance matrix C is ill-conditioned and numerical inversions 
are less stable. Due to broad covariance structure of the ma- 
trix C there are no completely independent information at 
different angular scales 6si, especially when there is no red- 
shift information. 



6.2 Use of redshift binning and third-order 
estimators 

In Figs. [T^17l we present our results from the Fisher anal- 
ysis described in sect. l6.l1 for various cosmological parame- 
ter pairs, combining the three angular scales 6s = 0.1', 10' 
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Figure 14. Results of Fisher matrix analysis for the parameter 
pair {nni,o"8}, combining the three angular scales 9s = 0.1', 10' 
and 1000'. Perfect knowledge of other parameters is assumed. 
For clarity only 3cr contours are displayed. Left panel: Fisher el- 
lipses obtained without redshift binning when we consider only 
the variance (M^p ) (estimator 1^.2 , dashed line) , only the third or- 
der cumulant (M^p)c (estimator Tis, long dashed line), and both 
cumulants for joint analysis (W2 and Wa, solid line). Right panel: 
Fisher analysis with two redshift bins. We consider the two bins 
separately for H2 (dashed and long dashed lines) and jointly for 
"Hi only (dot-dashed line) and for the joint pair and W3 (solid 
hne). 
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Figure 15. Same as Fig. I14l but for the {flm,^!!^ pair. 
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Figure 16. Same as Fig. I14l but for the {f2m,ns} pair. 



Figure 17. Same as Fig. 1141 but for the {as,ns} pair. 

and 1000' . Perfect knowledge of remaining parameters is as- 
sumed in each case and only 3a contours are plotted for clar- 
ity. We display both separate and joint analysis of second- 
order and third-order cumulants (estimators 0.2 and Hs) as 
well as the results obtained without redshift binning (left 
panels) and with two redshift bins (right panels). As ex- 
pected, we find that in all cases the 3a contours are much 
larger for the third-order cumulant {M^p)c than for the vari- 
ance {Mfp)c. Higher-order cumulants would be even more 
noisy. However, for the pairs {finufrg} and {fim, Wdc} the 
degeneracy directions (long axis of the ellipse) are signifi- 
cantly different so that combining Ti.2 and Ti.3 greatly im- 
proves the constraints on cosmology as compared with 7^2 
alone (see left panels) . Indeed, the long axis of the 7i2-ellipse 
can be fairly reduced by the intersection with the small axis 
of the 7^2-ellipse. This is most clearly seen in Fig. llSl for the 
{f2m,ii'do} pair. For the {^lni,ns} and {0-8,713} pairs where 
the 3a ellipses associated with Ti.2 and 7^3 are almost aligned 
the H2 contour is well within the TYa-ellipse so that adding 
the third-order cumulant does not tighten the constraints 
on cosmology. 

Next, we show in the right panels the effects of red- 
shift binning. For the {rim, era} pair the 7i2-ellipses have a 
very high eccentricity and are very thin so that the small 
change of orientation between the two redshift bins leads 
to a significant tightening of the intersection area which im- 
proves the constraints on these cosmological parameters. For 
other pairs of parameters the 7i2-ellipses associated with 
the lower redshift bin (1) are significantly broader than for 
the high redshift bin and do not improve the constraints 
on cosmology. Then, redshift binning is not really useful in 
these cases. The fact that the low redshift bin yields poorer 
constraints on cosmological parameters can be understood 
from the smaller line of sight which gives less room for weak- 
lensing effects. 

6.3 Comparison of Fisher analysis and 
likelihood 

We now compare the results of the Fisher matrix anal- 
ysis with the contour plots obtained through a likeli- 
hood function. We again combine the three angular scales 
9s = 0.1', 10' and 1000' and we only display the joint analysis 
of second-order and third-order cumulants. We first show in 
Fig. 1181 our results with no redshift binning. We can check 
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Figure 18. la and Scr Fisher ellipses compared with contours, 
combining the three angular scales ds = 0.1', 10' and 1000' with- 
out redshift binning. We only display the results corresponding 
to joint analysis of second-order and third-order cumulants. 
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Figure 20. Same as Fig. 1191 for other cosmological parameter 
pairs. 
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Figure 21. Same as Fig. 1191 for other cosmological parameter 
pairs including the typical source redshift zq- 



Figure 19. Same as Fig. 1181 but using tomography. The survey 
is split into two redshift bins and their constraints on cosmology 
are combined to yield the contours shown in the four panels for 
various parameter pairs. 

that the likelihood agrees well with the Fisher matrix 
analysis. Moreover, it happens that both Icr and So" con- 
tours are close to elliptic so that the Fisher matrix analysis 
appears to be sufficient, except for the broader 3ct contour 
of the {nmiWdc} pair which is large enough to see distor- 
tions from the elliptic shape (i.e. far from the fiducial model 
{Q,m ~ 0.3, Wde = —1} the deviations of cumulants {A'Ilp)c 
are no longer linear). Next, we show in Figs. 1191201 a compar- 
ison of Fisher matrix analysis with the likelihood func- 



tion using redshift binning. In agreement with sect. 16.21 we 
find that tomography is most efficient for the pair {ilm, os}- 
Both Fisher matrix and contour plots are consistent, as 
in Fig. 1181 Again, we find that the Scr areas are small enough 
to have elliptic shapes. 

6.4 Dependence on uncertainties on the mean 
redshift 

We now investigate the impact of any uncertainty on the 
source redshift distribution onto weak-lensing results. Thus, 
we show in Fig. 1211 Fisher matrix and contours which 
can be obtained for various cosmological parameters letting 
the typical source redshift zo defined in eg. lSlH free (in each 
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Figure 22. Marginalised 3-cr Fislier contours are displayed for 
various pairs of cosmological parameters. The inner solid contours 
correspond to the case when all hidden parameters are assumed 
to be known perfectly. Outermost long-dashed lines correspond 
to the case when all hidden cosmological parameters and median 
redshift zq are marginalised. The short-dashed line correspond to 
marginalisation of all hidden cosmological parameters except zq 
(i.e. a perfect knowledge of zq is assumed). In all cases two redshift 
bins are considered and information at the level of two-point and 
three-point are included for a SNAP class survey. Angular scales 
being considered are 8s = O.l', 10' and 1000' . 



panel other parameters are assumed to be known). Then, we 
can directly read in Fig. l21l bv how much any inaccuracy on 
the source redshift distribution can broaden the constraints 
on cosmology. Note that even when we do not impose any a 
priori on zq such weak-Iensing observations (assuming again 
that other parameters are given) are able to recover zq by 
themselves up to 20%. If zo is known to a better accuracy 
one can cut off the contours in Fig. 1211 Depending on the 
cosmological parameter of interest we can see that the un- 
certainty on Zo can multiply the error bar on the former by 
a factor of 2 (ris) up to 4 {Qm or as). Therefore, a good 
accuracy on the source redshifts can be quite rewarding. 



6.5 Marginalizing over unknown parameters 

Fig. 1221 shows the individual error ellipses for parameter 
pairs as various degree of marginalisation is used. For the 
case of two redshift bins and joint analysis of 2-point and 3- 
point analysis we plot the effect of having no prior knowledge 
of Zq or any other parameters on estimation error. Clearly 
accurate knowledge of power-spectrum and bi-spectrum evo- 
lution is essential for putting any constraints on spectral 
index ris or its run as . 

Results regarding marginalised errors are presented in 
units of f~i}J^ in tabled Although the entries represent an 
extrapolation of flat-sky results to spherical sky, it is un- 



likely that order of magnitudes are going to change with a 
more accurate all-sky approach. The second column shows 
Face for individual parameters where complete knowledge 
of all other parameters is assumed. The third column shows 
marginalised errors when other parameters and zo are un- 
known. Similarly the fourth column presents errors for addi- 
tional prior of knowing only zo perfectly. Last column corre- 
sponds to the determination of only Q^i, ag and dark energy 
equation of state w^c whereas all other parameters are as- 
sumed to be known. 

Finally, we present the cross-correlations Vij among es- 
timation errors for various parameter pairs in table |21 The 
signiflcant correlations show that it is difficult to measure 
simultaneously all parameters with a good accuracy, if there 
are no priors from other data sets. This is consistent with 
Fig. 1221 and with the comparison between the second and 
third columns of tabled 



7 DISCUSSION 

We have used weak lensing tomography in real space at the 
level of two-point (equivalently power spectrum) and three- 
point (equivalently bi-spectrum) to study how accurately 
the background dynamics and the nature of dark energy can 
be probed from future weak lensing surveys such as SNAP. 
It is well known that the weight functions and the growth 
rate dependencies are different for the second (M^p) and 



third-order (M^p) moments. These complimentary informa- 
tion help to reduce the level of degeneracies along with the 
tomographical information in certain speciflc choices of cos- 
mological parameters. 

In our formalism, cross-correlation among various an- 
gular scales at different redshifts can be very easily incor- 
porated in a natural way. Such a treatment is completely 
analytical and is an extension of previous study by Mun- 
shi & Coles (2002) which was later extended in Munshi & 
Valageas (2005). Assuming a "no-hole" approach it is pos- 
sible to directly include all contributions to covariance ma- 
trices which include cosmic variance at large angular scales, 
shot noise at small angular scales and mixed terms at inter- 
mediate scales. Taking advantage of a complete analytical 
description, we have check speciflc approximations to these 
covariances and their impact on error estimates of cosmolog- 
ical parameters. Over-enthusiastic simplifications of covari- 
ance structure of high-order estimators such as three-point 
estimators for (Mfp), which depend increasingly on accu- 
rate description of non-Gaussianities, are typically absent 
in harmonic domain based approach and can lead to erro- 
neous error-estimates. 

Assuming a WMAP-centric ACDM cosmology, we have 
considered several cosmological parameters such as Qui, fg, 
spectral index Us, running of spectral index Us, which de- 
termine the background dynamics of the universe, and the 
dark energy equation of state Wde- 

Tomography in some cases can help to reduce the level 
of degeneracy in {as, flm} or {w^e, Qm} however interest- 
ingly in a large number of other cases most of the infor- 
mation seems to be coming only from the higher redshift 
bin and sub-dividing the sources in different redshift bins 
does not seem to affect the estimation accuracy. Equiva- 
lently while inclusion of higher order moments does indeed 
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Table 1. Scatter in estimated parameters in units of f 



-1/2 
sky 



for individual and joint estimation. Angular scales involved are 0.1', 10' 



and 1000'. Values quoted within parenthesis are from two-point analysis whereas others are from joint analysis of two- and three-point 
statistics. Information from two redshift bins is considered. Different columns show the cases where we know all other parameters (2nd 
column), we marginalize over other cosmological parameters Xj and redshift zq (3rd column), we marginalize over other parameters Xj 
only (4th column) and only the three parameters {Oni, o"8> '"'del ^^e unknown. 
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Table 2. Cross-correlation Tij among estimation errors for various parameters. Values quoted within parenthesis are from two-point 
analysis whereas others are from joint analysis of two-point and three-point statistics. Two redshift bins are considered. Angular scales 
involved are 0.1', 10' and 1000'. 







o"8 


Wdc 


Us 


as zq 




-t-1.00 










0-8 


-0.88( -0.96) 


+1.00 










-l-0.96(-f0.60) 


-0.92(-0.37) 


+1.00 






ris 


-0.87(-0.97) 


-f0.56(-f0.90) 


-0.78(-0.63) 


+1.00 




as 


-l-0.75(+0.59) 


-0.37(-0.53) 


+0.64(+0.29) 


-0.97(-0.75) 


+1.00 


ZQ 


-l-0.92(+0.99) 


-0.94(-0.98) 


+0.97(+0.53) 


-0.73(-0.95) 


+0.58(+0.55) +1.00 



help to break the parameter degenerecy for some parame- 
ter pairs, the other combinations are largely unaffected by 
inclusion of third order information. Clearly both tomogra- 
phy and non-Gaussianity information are more useful when 
a joint estimation of several parameters is performed. This 
is particularly true for spectral index Us and its running as 
which enter the modelling of bi-spectrum only through the 
description of power spectrum. Previous studies of tomo- 
grpahy made use of Monte-Carlo simulations of correlated 
field of views (see e.g. Simon et al. 2004) across redshift bins. 
Our analytical results which do not rely on numerical simula- 
tions validate and generalise similar studies by including de- 
tailed descriptions of non-Gaussianities and cross-correlatios 
among various redshift bins. 

Generalised third order moments were considered re- 
cently by Kilbinger & Schneider (2005) which seem to be 
improving the situation by adding extra information at the 
level of third order. Interestingly however the highly cor- 
related nature of such additional input suggests that one 
would only be restricted to a limited range of angular scales 
for construction of such a generalised quantities. A more 
general discussion of such generalised statistics which are 
also named as cumulant correlators in the literature can 
be found in Munshi & Valageas (2005). In our calculation 
here we have ignored the primordial non-Gaussianity which 
can potentially be studied at low redshift given the large 
sky-coverage of future weak lensing surveys. Recent results 
by Takada & Jain (2005) showed lensing fields do retain 
some information regarding primordial non-Gaussianity al- 
though it is expected that most of the information regarding 
baryon oscillations and primordial non-Gaussianities will be 



erased by large scale non-linearities generated by gravita- 
tional clustering at a later epoch. Extension of our results 
in such directions to take into account additional cosmolog- 
ical parameters, with priors from external data sets such as 
Type-IA supernova or CMB observations, will be presented 
elsewhere. Additional parameters related to time evolution 
of dark energy equation of state or neutrino mass will also be 
included. Despite large sky-coverage even with future space 
based missions such as SNAP, error of estimation in vari- 
ous cosmological parameters remain inherently degenerate 
due to very nature of observables in weak lensing surveys. 
This is particularly underlined by the very high values of 
degradation factor (Faa)~^''^/(F~„)^''^ for various levels of 
marginalisation as presented in table-0 In table-|5|scaled ver- 
sion of inverse Fisher matrix rafs which is also the cross cor- 
relation coefiicient r^^ = (F-i)^^/{(F-^)cc(F-i)^^}i''2 of 
estimation error of different parameters is plotted for cer- 
tain choices of the parameters as indicated. Our findings 
are in agreement with a recent study by Kilbinger & Mun- 
shi (2005) where they reached a similar conclusion based 
on more elaborate generalised eigenmode analysis of Fisher 
matrices. 

Knowledge of redshift distribution of source galaxies 
from external observations can help parameter estimation 
especially for smaller surveys. Even for larger surveys such 
as SNAP we found that order of magnitude improvement in 
accuracy is possible for almost all parameters we have stud- 
ied. Large scale spectroscopy of a fair sample of faint source 
galaxies may not be possible, alternatively photometric red- 
shift determination of subsample of observed galaxies can 
still improve the level of accuracy. It is likely that future ob- 
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servations with near complete sky coverage fsky ~ 1 span- 
ning large redshift coverage can deal with much more de- 
tailed description of the source redshift distributions which 
we have only parametrised by one parameter zo- 

Several authors have recently studied joint constraints 
on various cosmological parameters such as {ns,<Js} and 
{qs , Ua } using Lyman-alpha forest data jointly with 1-year 
observations from WMAP (Spergel et al. 2003). Clearly the 
future space based SNAP class experiments will help to pro- 
vide competitive constraints when used jointly with external 
data set. 

Bi-spectrum of density fluctuations carries complemen- 
tary information regarding background geometry and dy- 
namics of the universe. Although current surveys are at the 
detection limit (see Bernardeu 2002) of measuring significant 
non-Gaussianity, it is expected that future space based ob- 
servations with much greater sky coverage will supplement 
power spectrum information with accurate measurement of 
bi-spectrum. As we have shown here, this will require a very 
accurate description of evolution of bi-spectrum. Associat- 
ing unknown parameters such as /2 and fs with low order 
description of dark matter clustering and trying to estimate 
them directly from the data degrades the estimation accu- 
racy of other cosmological parameters considerably. 

Finally, in our studies we have compared Fisher analy- 
sis with full grid based calculations. Importantly Fisher 
analysis should only be seen as a way to estimate minimum 
variance errors for unbiased estimators and to check degener- 
acy directions. It does not propose to reveal the detailed be- 
haviour of the confidence plane far from the flducial models. 
This is particularly true for very asymmetric constraints e.g. 
in (fimjCg) plane. Nevertheless the various parameter pairs 
that we have studied show a reasonable match in confidence 
level at la and even to 3a levels. An alternative approach 
which uses entropy functional to map out the likelihood dis- 
tribution in the parameter space was considered in Taylor 
& Watts (2002). A grid based calculation of full becomes 
prohibitively costly as the dimension of the parameter space 
increase. A more effective sampling of the parameter space 
can be achieved by Monte Carlo based Markov-Chain algo- 
rithms. A full implementation will be presented elsewhere. 

In our present analysis we have focused on survey pa- 
rameters similar to proposed SANP-class surveys (JDEM). 
It is not difficult to extend it to other survey speciflcations. 
With future well optimised surveys with large sky cover- 
age, CFHTLS (172 deg^), SNAP (300 deg^), VISTA (10000 
deg'^), Pan-starr (31000 deg^) and recent data analysis tech- 
niques, weak lensing will be a very useful tool to probe not 
only standard cosmological parameters, such as fim and as, 
but also the dark energy equation of state parameters such 
as Wdc- However to achieve this goal it is essential to have a 
tight handle on systematics which should at least stay com- 
parable to statistical errors. 
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APPENDIX A: OTHER SURVEY 
CHARACTERISTICS 

For completeness, we display in this appendix the noise to 
signal ratios we obtain for different survey properties. Thus, 
we display in Fig. lAll our results for the SNAP survey as in 
Fig. El (solid lines) as well as for the cases when we divide 
by a factor two its area (dashed line), its galaxy number 
density (dotted lines) or its depth (dot-dashed lines). When 
the depth of the survey is varied the dashed curves labeled 
/2 and /s change slightly (as the redshift distribution of 
sources whence of the density fluctuations probed by the 
survey is varied) but since the modiflcation is very small (by 
definition the curves show the same plateau at small scales) 
we do not add them in Fig. lAll for clarity. We can check that 
a smaller galaxy number density increases the noise at small 
scales (below 2') which is dominated by shot noise but does 
not change the signal to noise ratio at large scales which 
is governed by the cosmic variance. By contrast, a smaller 
area leads as expected to a smaller signal to noise ratio at all 
scales. Finally, a smaller depth decreases the signal to noise 
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Figure Al. Inverse of the signal to noise ratios of the estimators 
'H2 and Tis of the second-order and third-order cumulants of the 
aperture-mass for various survey properties. The solid lines corre- 
spond to the SNAP survey as in Fig.|2l the dashed lines to a sim- 
ilar survey with an area twice smaller (A = 150deg^), the dotted 
lines to a galaxy number density twice smaller (rig = 50arcmin~'^) 
and the dot-dashed lines to a survey with a depth which is twice 
smaller (zq = 0.57 and Zmax = 1.5). There is no redshift binning 
and we use the full covariance Ca. 
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Figure A2. Inverse of the signal to noise ratios of the estimators 
W2 and Tis of the second-order and third-order cumulants of the 
aperture-mass as in Fig. |21but for a typical ground-based survey. 
There is no redshift binning. The solid lines correspond to the 
full covariance C while the dotted line show the results we obtain 
when we only use its Gaussian part . 



ratio at small scales as the amplitude of gravitational lensing 
effects is weaker because of the smaller line of sight while 
the shot noise remains the same. At large scales the signal 
to noise ratio for 7^2 is not diminished while it improves 
somewhat for Ti.3 (at fixed galaxy number density) as the 
survey focuses more on the non-linear regime of gravitational 
clustering where non-Gaussianities are larger. 

Next, we show in Fig. IA2l the noise to signal ratios ob- 
tained for a typical ground-based survey with a galaxy dis- 
tribution n{za) (X z'^e~^^'^°'^^ , an area A = 10 deg^ and a 
galaxy immber density Ug — 20 arcmin"^ (e.g., van Waer- 
beke et al. (2001)). Of course, we can check that both the 



shot-noise and the cosmic variance are significantly larger 
than for a space-based mission such as SNAP. This implies 
that the range of scales which can be used to extract cos- 
mological information is quite smaller. On the other hand, 
the theoretical uncertainty on the non-linear regime of grav- 
itational clustering (denoted by the dashed curves) is only 
of the same order as other sources of noise around 5' and is 
relatively smaller at other scales, so that it is not really the 
limiting factor for such surveys. 

A generalised cigcii mode analysis (Locvc 1948, 
Karhunen 1947) which can address the issue of optimisa- 
tion of survey design [zo, f sky, rig) for recovery of a given set 
of cosmological parameters in line with Kilbingcr & Munshi 
(2005) but including tomography and non-Gaussianity will 
be presented elsewhere. 
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